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Magnetization of the frustrated S — 1/2 chain compound LiCuV04, focusing on high magnetic 
field phases, is reported. Besides a spin-flop transition and the transition from a planar spiral to 
a spin modulated structure observed recently, an additional transition was observed just below the 
saturation field. This newly observed magnetic phase is considered as a spin nematic phase, which 
was predicted theoretically but was not observed experimentally. The critical fields of this phase 
and its dM/dH curve are in good agreement with calculations performed in a microscopic model 
(M. E. Zhitomirsky and H. Tsunetsugu, preprint, arXiv: 1003.4096 ^2). 
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Unconventional magnetic orders and phases in frus- 
trated quantum spin chains are attractive issues, because 
they appear under a fine balance of the exchange inter- 
actions and are sometimes caused by much weaker inter- 
actions or fluct nations ji"— 

A kind of frustration in quasi one-dimensional (ID) 
magnets is provided by competing interactions when the 
intra-chain nearest neighbor (NN) exchange is ferromag- 
netic and the next-nearest neighbor (NNN) exchange is 
antiferromagnetic. This type of exchange results in the 
formation of bound magnon pairs in a high magnetic field 
near the saturation point. ^■^■*' The condensation of bound 
magnon pairs could provide a nematic spin ordering. Nu- 
merical investigations of frustrated chain magnets with 
different models^^— have, indeed, predicted a nematic or- 
dering of transverse spin components at zero temperature 
in a magnetic field range close to the saturation field. The 
spin nematic phase was introduced phenomenologically^*' 
as a magnetic state with zero ordered spin components 
and nonzero spin correlations emerging via a second or- 
der phase transition. This spin ordering is analogous to 
the well known ordering of molecules in nematic phases 
of liquid crystals. ^^ Besides the frustrated spin S = 1/2 
chains, S ^ 1 systems with biquadratic exchange were 
proposed as potential nematics.^^ Nevertheless, the real 
spin nematics have not been found during a long period. 

Recently a new analytical approach for the description 
of the high field magnetic phase of frustrated chains with 
a spin- nematic ordering was proposed. ^'^ The parameters 
of this curious phase were calculated from the represen- 
tation of the condensate wave function in terms of the 
coherent bosonic states. A S = 1/2 quasi ID magnet 
with the competing interactions of the described type, 
LiCuV04, was analyzed as a potential spin nematic. For 
LiCuV04, the nematic phase was predicted to exist in 
the field interval of about 3 T below the saturation field. 
This work stimulated our high field magnetization mea- 



surements on LiCuV04 single crystals. 

LiCuV04 crystallizes in an inverse spinel structure 
AB204 with an orthorhombic distortion. Copper ions 
share the octahedral sites being arranged in chains di- 
rected along the 6-axis. These chains are separated by 
the nonmagnetic ions of Li, V, O. The elementary cell 
contains four ions Cn^~^{S = 1/2) (see Fig. [J). Neutron- 
diffraction experiments^ indicate that in the ordered 
phase at r < T/v = 2.3 K and H = an incommen- 
surate planar spiral spin structure occurs with the wave 
vector kic — (0.468) • 27r/5 directed along the 6-axis. The 
Curie- Weiss temperature is Qcw = — 15 K. The ordered 
magnetic moments of Cu^+ ions lie within the a6-planes 
and amount to 0.31/ib at T=1.6 K. The strong reduction 
of the ordered component and a small value of T^ /Qcw 
indicate a significant influence of quantum fluctuations 
and frustration on the ground state. 

Inelastic neutron-scattering experimentsi^ confirmed 
that the exchange interactions Ji_2 within chains are 
dominant, and the NN and NNN intrachain exchange 
interactions are ferromagnetic and antiferromagnetic, re- 
spectively. It was shown that the incommensurate struc- 
ture is due to these competing interactions. In addi- 
tion, it was found that there is a ferromagnetic exchange 
between neighboring chains along the a-axis (J3), and 
an antiferromagnetic exchange link between different ab- 
planes ( J4), which is almost five times smaller than J3. 

In a magnetic field H, consecutive magnetic phase 
transitions were observed f^^^-i^ Below 12 T, the magnetic 
phases were studied also by means of ESR and NMR tech- 
niquesJ^ For H within the a6-plane, the first transition 
takes place at iiqHci ~ 2.5 T, being interpreted as a spin- 
flop reorientation of the spin spiral plane perpendicular 
to HJ^ At noHc2 « 7 T, there is a transition to a phase 
with the coUinear spin-modulated structure. The spins 
in this phase are directed along the magnetic field. This 
ordering occurs only within separate a&-planes, and ad- 




FIG. 1: (Colour online) Scheme of Cu^^ moments in the crys- 
tal structure of LiCuV04. Copper ions are marked by circles. 
Arrows constitute a helical spin structure at H — below Tn 
(Ref. Il4l'). Ji, J2, and J3 are main exchange integralsl^. The 
exchange bonds Ji, J2, J3 lie in shadowed planes, which are 
parallel to the ab-planes. 



jacent planes are not correlated . ^^1^^ A week anisotropy 
of Hc2 reveals an exchange nature of this transition. A 
possible mechanism of this transition was proposed in 
Refs. 13 and |^. Further, the saturation of the magne- 
tization was ascribed to dM/dH anomalies observed at 
fJ-oH~40 T and 45 T for the orientations ^ || c and H \\ a, 
respectively] ^^'^^ However, the true saturation, i.e. the 
transition to constant M, was not achieved there. 

Single crystals of LiCuV04 of several cubic millime- 
ters were prepared as described in Ref. 1191 . The com- 
position of samples is characterized by the relations: 
Li/V=0.96±0.05, Li/Cu=0.95±0.04, Cu/V=0.99±0.01. 
Thus, the average composition is Lio.97CuV04, which 
means a low concentration of vacancies in the Li- 
sublattice. This type of crystals has the unit cell param- 
eter c=8.745 A. The analysis of the quality of samples 
is discussed in Ref. |20. The samples used in the present 
work are from the butch 1 of Ref. |20. Magnetization 
curves in magnetic fields of up to 7 T were measured 
with a commercial SQUID magnetometer (Quantum De- 
sign MPMS-XL7). High field magnetization data were 
taken in magnetic fields of up to 55 T using a pulsed 
magnet at KYOKUGEN in Osaka University^i 

Experimental dM/dH ioi H \\ c curve and its field 
integral M{H) for H \\ a, 6, c are shown in Fig. 2. Each 
dM/dH curve demonstrates two sharp peaks at the fields 
Hc2, Hc3 and a step at Hsat- A strong inflection is seen 
on the M{H) curve below the field Hcs- The lowest tran- 
sition, marked as Hc2, occurs at iJ « 7 T for all three 
field directions. This anomaly corresponds to the tran- 
sition from a spiral to a collinear spin modulated phase 
described above. The field derivative of magnetization 
above Hc2 is a bit smaller than that in the low field phase. 
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FIG. 2: (Color online. )Upper panel: dM/dH curve measured 
with a pulse field technique. Insert: dM/dH near Hsat- The 
solid line presents the fit by a step-like function. The dashed 
line marks zero level. Lower panel: M{H) curves at T =1.3 
K obtained by integration for H \\ a,b,c. The solid lines 
correspond to the magnetization for field descending process. 
The straight dashed lines present linear fits in the range Hc3 < 
H < Hsat and extrapolations of these fits to zero field. 



This difference equals 5±1% for i/ || c, 3.5±1% for iJ || a 
and 2.5± 1% for H \\ b. At this critical field, a jump of 
magnetization for 14±2% occurs. 

Our intriguing finding in the magnetization curve is a 
"fine structure" near the saturation process, which in- 
cludes anomalies at Hc3 and H^at- The field Hc3 is 
marked by the maximum of dM/dH. Right under Hcs 
(Fig. 2), the dM/dH increases abruptly more than by a 
factor of twenty. The magnetic moment at Hc3 reaches 
95% of the saturation value. Above Hcs, the value of 
dM/dH drops quickly to that in the middle field range. 
Then, it goes down to zero at the saturation field H^at- 
The value of Hgat is determined by the fit of dM/dH 
curve in the interval above Hcs by a Fermi-like step func- 
tion: dM/dH = x/{exp{{H - Hsat) /AH) + 1) where 
Hsat, the derivative dM/dH = x, and the width of the 
saturation transition AH are fitting parameters. The ob- 
tained value of Hsat does not depend on the boundaries 
of the fitting range within 1%. The value of x is ob- 
tained with the accuracy of 20%, AH is less then 0.01 T. 
The magnetic moment in the saturated phase for H\\c is 
M — 1.15 ± 0.05 /is/Cu^"^, consistent with the expected 
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FIG. 3: (Colour online.) dM/dH, M{H) curves measured 
with a pulse field technique for H\\c. The dashed and solid 
lines correspond to the magnetization curves for field ascend- 
ing and descending processes, respectively. 



1.15 ^B for gc—2.3. The fields Hcs and Hgat are well sep- 
arated by an interval of 3.9 ±0.1 T for H \\ c, 4.9±0.2 T 
for H II a, b. 

The temperature evolution of the M{H) curve for 
i/ II c is shown in Fig. 3. The small kink of M{H) at 
Hsat is distinguishable below 2 K. Above 2 K, the phase 
transition at i/c3 is broadened and the anomaly at Hsat 
is masked. At 1.3 K, the hysteresis is practically negligi- 
ble, while at higher temperatures it becomes larger and 
is probably caused by the magnetocaloric effect. 

The results of low field measurements are depicted in 
Fig. 4. The dM/dH curves for H \\ a and H \\ b have 
peaks at Hci=2.5 T and 3.1 T, respectively, correspond- 
ing to the spiral-plane flop transition. The critical fields 
are listed in Table 1. The proposed magnetic phases at 
T < 2 K are presented in the lower panel of Fig. 4. For 
-ff II c, the low field transition is absent, because the 
spiral-plane flop does not take place in this case. The val- 
ues of Hci for a and b directions arc consistent with those 
obtained from ESR measurements.^^ The Hc2 for H \\c 
is close to that observed in the dM/dH curve in Ref. [l3- 
The field Hc3 corresponds well to the peak of dM/dH 
observed in Ref. llSl where this anomaly was ascribed to 
the saturation field. Nevertheless, the right wing of this 
dM/dH peak at the high field range of their work ap- 



peared much more diffuse than that in the present obser- 
vation. The interval where the dM/dH decreases twice is 
for this experiment about 2 T while for our observation 
this interval is 0.25 T. A large width of the right wing 
of the dM/dH peak masks the details of saturation in 
the experiments of Refs. [Ij and [la. Our measurements 
make clear that the field at the large peak of dM/dH 
is not the saturation field Hgat as supposed before, but 
marks a transition into a new high field phase. 
For the analysis, we use the model Hamiltonian: 

'H = 2_J[Jl{^i,j ' Si,j+l) + Si^j ■ J I ■ Si^j+i] (1) 

+ / \J2\Si.j ■ Sij^2) + Si,j ■ J 2 ■ Sij^2] 

+ ^I^bH g ■ S^.j, 

where gaa = 2.070, gbb = 2.095, and gcc = 2.313 are 
the diagonal components of 5-tcnsor, Ji/Zcb — —18.5 K, 
Ja/fcB = 44 K, Js/fcB = -4.3 K and Jf.JkB = 
0.16 K, Ji^bh/fcB = -0.02 K, and Jf^^/ZcB = -1-75 K, 
"^Icc/^B — —0.2 K are the diagonal components of J° 
and J2 -tensors. The remaining components of tensors 
are proposed to be small. The summation is over ex- 
change bonds marked in Fig. 1. The components of g 
and J° tensors are obtained from ESR experiments^^, 
"^1,2,3 ~ from the neutron scattering studyiS, and J2cc 
- from the gap of the low-frequency antiferromagnetic 
resonance mode^-. The saturation fields H^^^''^ may 
be evaluated roughly in the classical ID approximation 
as iJJ,t=4J25(l + Ji/AJ2f/g,^iB ii=a,b, and c)^^^: 
tioH^,t=50^ T, Moi?Lt=50 T, and Moi?.^at=45 T. These 
values are in good agreement with our experiment. 







TABLE I: 








i/ci (1.9 K) 


Hc2 (1.3 K) 


^c3 (1.3 K) 


Hsat (1.3 K) 


H\\a 
H\\b 
H\\c 


2.5 T 
3.1 T 


7.6±0.2 T 
7.6±0.2 T 
6.7±0.2 T 


47.9±0.2 T 
47.2±0.2 T 
40.5±0.2 T 


52.8±0.3 T 
52.1±0.3 T 
44.4±0.3 T 


Theory" 
(T = 0) 


- 


- 


44.5 T 


47.1 T 



Our observations clearly indicate that the field Hc3 is 
different from the saturation field. The field and mag- 
netization intervals, corresponding to the newly found 
magnetic phase, can be compared with theoretical calcu- 
lations performed in different models for a nematic phase. 
A numerical investigation of a frustrated chain was car- 
ried out^ using a model with anisotropic exchange in- 
teractions Jfcc, J2cc: at the condition Jicc/Ji=J2cc/J2- 
Three phases were predicted for a system with exchange 
parameters Ji, J2, and Jj^^^ of LiCuV04. In low fields, 
there is a quadrupolar spin density wave phase, having 
the same nearest spin correlation as the spin modulated 
phase observed in the experiment. A quadrupolar spin 
nematic phase should occur in a range below H^at ■ Then, 
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FIG. 4: (Colour online) Upper panel: the low field magne- 
tization curves M{H) (left) and their field derivatives (right) 
for three field directions and at T^l.9 K. Lower panel: the 
scheme of expected magnetic phases in LiCuV04 at T < 2 K. 



a saturated phase appears above Hsat- This work pre- 
dicts a nematic phase in the interval of magnetization 
between Mgat — AM and Mgat ■ The value of AM is pre- 
dicted to be 20-30 % of Mgat ■ The experimental value of 
AM in LiCuV04, however, is 5% of Mgat- In low fields, 
a planar spiral is observed for LiCuV04 instead of a spin 
density wave predicted by the model. Note that the ex- 
perimental M{H) curves differ essentially from those ob- 
tained numerically for a ID model^. For ID model, the 
fluctuation divergence of dM/dH takes place at Hsat, 
whereas in the experiment the divergence is found at 
Hcz < Hsat- Probably, this difference is caused by the 
2D character of exchange links in LiCuV04. 



A more adequate 2D model with exchange interactions 
Ji, J2, J3 was considered analytically in Ref.[lJ. A phase 
transition from a spiral umbrella-like phase to a spin- 
nematic one was predicted at a critical field He- The 
field range of the nematic phase should reach Hsat- For 
exchange parameters of LiCuV04, the saturation field is 
predicted to be 47.1 T, and the width of field range of 
nematic phase Hsat — He— 2.6 T, which is not far from 
the experimental value of Hsat-Hcz (3.9 - 4.9 T). This 
theory presumes an isotropic ^-factor of 2.0, which dif- 
fers from the real values of ga,b,c- Besides, the exchange 
parameters Ji,2,3 are known to have the accuracy of 12, 
1.5 and 70%, correspondingly. This may be a reason of a 
considerable uncertainty in the estimation of He- More- 
over, we compare the value of dM/dH in the high field 
phase with the theoretica.1 prediction for the spin nematic 
phase. According to Ref . llSl. dM/dH in the nematic state 
is approximately a half of the slope of the classical mag- 
netization curve, which is Msat/ Hsat- As one can see in 
Fig. 2, the value of dM/dH in the high field phase corre- 
sponds well to ^Msat/Hsat indeed. 

In conclusion, a new magnetic phase near the satu- 
ration field is observed in a quasi- ID S — 1/2 mag- 
net with competing ferromagnetic NN and antiferromag- 
netic NNN exchange interactions LiCuV04. By com- 
parison with theoretical studies of Refs. Il]|7l-l9lfl3l. this 
new phase must be a spin nematic one. This proposition 
is supported by reasonable agreement between the criti- 
cal fields and dM/dH calculated analytically^"^ in the 2D 
model and experimental ones. 
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